Capacity of nonlinear bosonic systems 
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We analyze the role of nonlinear Hamiltonians in bosonic channels. We show that the information 
capacity as a function of the channel energy is increased with respect to the corresponding linear 
case, although only when the energy used for driving the nonlinearity is not considered as part of 
the energetic cost and when dispersive effects are negligible. 
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Noninteracting massless bosonic systems have been the 
object of extensive analysis ^Ig^]. Their maximum capac- 
ity in transmitting information was derived for the noise- 
less case both in the narrow-band regime (where only few 
frequency modes are employed) and in the broad-band 
regime. However, it is still an open question whether 
nonlinearities in the system may increase these bounds: 
linearity appears to be the most important assumption in 
all previous derivations 0- Up to now nonlinear effects 
have been used in fiber optics communications to over- 
come practical limitations, such as using solitons to beat 
dispersion or traveling wave amplifiers to beat loss [sj. 
The approach adopted in this paper is fundamentally dif- 
ferent from these and from others where nonlinearities 
and squeezing are employed at the coding stage when 
using linear channels jjf. We follow the cue of a recent 
proposal where interactions were exploited in increas- 
ing the capacity of a qubit-chain communication line, fn 
the case of linear bosonic systems, the information stor- 
age capacity of a signal divided by the time it takes for it 
to propagate through the medium gives the transmission 
capacity of the channel. In the presence of nonlineari- 
ties, dispersion can affect the propagation of the signal 
complicating the analysis, but an increase in the capac- 
ity can be shown, at least when the dispersive effects are 
negligible. A complete analysis of dispersion in nonlinear 
materials is impossible at this stage, since the quantiza- 
tion of these systems has been solved only perturbatively. 
The basic idea behind the enhancement we find is that 
the modification of the system spectrum due to nonlinear 
Hamiltonians may allow one to better employ the avail- 
able energy in storing the information: we will present 
some examples that exhibit such effect. Excluding the 
down conversion channel (a model sufficiently accurate 
to include the propagation issue -see Sec. Iff C|) . all these 
examples are highly idealized systems but are still indica- 
tive of the possible nonlinearity-induced enhancements 
in the communication rates. An important caveat is in 
order. In the physical implementations that we have an- 
alyzed, there is no capacity enhancement if we include 
in the energy balance also the energy required to create 
the nonlinear Hamiltonians. This is a general charac- 
teristic of any system: if one considers the possibility of 
employing all the available degrees of freedom to encode 



information, then one cannot do better than the bound 
obtained in the noninteracting case However, the 

enhancement discussed here is not to be underestimated 
since in most situations many degrees of freedom are not 
usable to encode information, but can still be employed 
to augment the capacity of other degrees of freedom. A 
typical example is when the sender is not able to modu- 
late the signals sufficiently fast to employ the full band- 
width supported by the channel: an external pumping 
(such as the one involved in the parametric down conver- 
sion case) may allow to increase the energy devoted to 
the transmission modes. 

We start by describing a general procedure to evalu- 
ate the capacity of a system and we apply it to a lin- 
ear bosonic system to reobtain some known results (see 
Sec. nj. Such a procedure is instructive since it empha- 
sizes the role of the system spectrum in the capacity cal- 
culation. We then analyze a collection of examples of 
nonlinear bosonic systems in the narrow-band and wide- 
band regimes (see Sec. 1^): for each case we describe 
the capacity enhancements over the corresponding linear 
systems. General considerations on the energy balance 
in the information storage and on the information prop- 
agation conclude the paper (see Sec. IIII|) . 



I. INFORMATION CAPACITY 

The information capacity of a noiseless channel is de- 
fined as the maximum number of bits that can be reliably 
sent per channel use. From the Holevo bound (proven in 
Ref. [l| for the infinite dimensional case), we know that 
it is given by the maximum of the Von Neumann entropy 
S{g) = —TT[g\og2 g] over all the possible input states 
g of the channel. In our case, since g is the state of a 
massless bosonic field, the associated Hilbert space is in- 
finite dimensional and the maximum entropy is infinite. 
However, for all realistic scenarios a cut-off must be intro- 
duced by constraining the energy required in the storage 
or in the transmission, e.g. requiring the entropy S{g) to 
be maximized only over those states that have an average 
energy E, i.e. 



E TT[gH] 
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where H is the system Hamiltonian. The constrained 
maximization of S{p) can be solved by standard varia- 
tional methods (see for examples), which entail the 
solution of 



6 S{g) 
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where A and A' are Lagrange multipliers that take into 
account the energy constraint and the normalization 
constraint Tr[g] — 1, and where the In 2 factor is in- 
troduced so that all subsequent calculations can be per- 
formed using natural logarithms. Eq. Q is solved by the 
density matrix g = exp[— AiJ]/Z(A) where 



Z{X) = Tr[e-^^] 



(3) 



is the partition function of the system and A is deter- 
mined from the constraint Q by solving the equation 



d_ 
dX 



InZ(A) . 



(4) 



The corresponding capacity is thus given by 



C^S [exp(- Aff ) /Z(A)] = [XE + In Z(A)] / In 2 , (5) 

which means that we can evaluate the system capacity 
only from its partition function Z{X). 

In general an explicit expression for Z(X) is difhcult 
to derive, but it proves quite simple for noninteracting 
bosonic systems, such as the free modes of the electro- 
magnetic field. In fact, in this case the Hamiltonian is 
given by 



H 



alak 



(6) 



which is valid in the limit HSuj ^ E. When applied 
to communication channels, this last equation is usually 
expressed in terms of the rate R (bits transmitted per 
unit time) and power P (energy transmitted per unit 
time) as 



R 



(10) 



by identifying the transmission time with 2Tr /Slo. In the 
next section we analyze how the capacities Cnb and Cu,b 
are modified by introducing nonlinear terms in the sys- 
tem Hamiltonian. 



II. NONLINEAR HAMILTONIANS 

Nonlinear terms in the Hamiltonian of the electro- 
magnetic field derive from the interactions between the 
photons and the medium in which they propagate. In 
this section we will employ the techniques described 
above to derive the narrow-band and wide-band ca- 
pacities when quadratic nonlinearities are present. In 
Sees. Ill Al IIIBI and III CI we discuss parametric down- 
conversion type Hamiltonians in the narrow-band and 
broad-band regimes. In Sees. Ill Dl and UTeI we discuss a 
mode swapping interaction. All these nonlinearities arise 
in real- world systems from x^^^ type couplings, when one 
of the three fields involved in these kinds of interactions 
is a strong pump field that can be considered classical 
[To| . For all the cases analyzed we present the capac- 
ity enhancement over the corresponding non interacting 
Hamiltonian. 



where ujk is the frequency of the fcth mode and the 
mode operators satisfy the usual commutation rela- 
tions [afc,aj^,] = Skk'- Hence, the partition function is 



z{x) = T] 



oo 

k rifc=0 
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From Eqs. (0 and 10 it is clear that the capacity C 
will be ultimately determined by the spectrum ujk of the 
system. In particular, in the narrow-band case (where 
only a mode of frequency lo is employed) Eq. (|3J) gives 0, 



(8) 



where g(x) = (1 + a;)log2(l + x) — xlogjX for x ^ 
and g{0) — 0. On the other hand, in the homogeneous 
wide-band case (where an infinite collection of equispaced 
frequencies uik = kSuj is employed for k e N) Eq. (0) 

gives Hasg 



wh 



TT 
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A. Squeezing Hamiltonian 

Consider the single mode described by the Hamiltonian 

H ^ huja^a + [(at)2 + a^] /2 + hn{^) , (11) 

where ^ is the squeezing parameter. We employ |^| < uj 
to avoid Hamiltonians which are unbounded from below. 
In Eq. the frequency O(^) = ^(w - y^oj^ - ^^) has 
been introduced so that the ground state of the system is 
null: with this choice, the average energy E is the energy 
associated to the mode a in the nonlinear medium. By 
applying the canonical transformation 



= A cosh 6 - A'^ sinh ( 
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(12) 



the Hamiltonian H is transformed to the free field form 
hy^uj^^^^A^ A, so that the derivation of the previous 
section can be employed to calculate the capacity. It is 
thus immediate [see Eq. ((HJ] to find the capacity of this 
system as 



C^g 



E 



(13) 
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This quantity measures the amount of information that 
can be encoded when E is the total energy associated 
with the mode a, and a nonhnear squeezing-generating 
term is present in the Hamiltonian: this result is quite 
different from the one obtained by using squeezed states 
as inputs to a linear system (see for example 0). The 
capacity C of Eq. (|13|l is higher than the capacity Cnb of 
the linear case ^ = since g{x) is an increasing function 
(see Fig.nj. The reason behind this enhancement is that 
the nonlinearity reduces the effective frequencies of the 
modes (from uj to y'uj^ — ^) , so that more energy levels 
can now be populated with the same energy. 



C 




FIG. 1: Capacity increase of the squeezing Hamiltonian of 
Eq. IIH as a function of the energy parameter E/fiu and 
of the squeezing ratio ^/uJ. The increase is evident from 
the positivity of C — Cnb, where C is given in Eq. 1131 and 
Cnb is the free space narrow-band capacity of Eq. No- 
tice that for high energy, the increase tends asymptotically 
to -log2 [1-^V^'] /2. 



B. Two mode parametric down conversion 

Consider the two interacting modes a and b evolved by 
the Hamiltonian 

H = hLu{a^a + b^b) + hC{a''b^ + ab) + , (14) 

where |^| < is the coupling constant and f2(^) = lo — 
\J llP' — has again been introduced to ensure that the 
energy ground state is null. [This Hamiltonian, as the 
previous one, possesses the structure of the algebra of the 
SU(1,1) group, in the Holstein-Primakoff realization] . We 
can again morph the Hamiltonian to the free field form 
using the two-field canonical transformation 



A — a cosh d 
B = a^ sinh 6* 



V b cosh d 



= - In 

4 



,(15) 



which transforms the Hamiltonian 
h^fZ]^^^{A^A- 



((Tl|l to the form 
B^B). The capacity in this case is 



C^2g 



E 



(16) 



which is higher than the capacity of two independent 
single-mode bosonic systems, given by Eq. I|16|l for ^ = 
0. [The factors 2 in p6|l derive from the presence of 
the two modes A and B.] The enhancement is again a 
consequence of the fact that the nonlinearity reduces the 
effective frequency of the two modes: lu — > -^/o;^ — 



C. Broadband parametric down-conversion 

Here we apply the above results to the case of two 
wide-band modes (signal and idler) coupled by a para- 
metric down conversion interaction. The interaction is 
mediated by a nonlinear crystal with second order sus- 
ceptibility x^"^"^ pumped with an intense coherent field 
of amplitude £p at frequency ujp. Assuming undepleted 
pumping and perfect phase matching, the Hamiltonian 
at the first order in the interaction is given by 

H = ^2 ^f^ko-lak + h{LUp - LUk)b\bk 

k 

+Kk{albl + akbk) + , (17) 

where Ofe and bk are the mode operators of the downcon- 
verted modes of frequency LUk and LUp — LUk respectively. 
Their interaction is described by the coupling parameter 
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ceo UaiuJk) nb{ujp - UJk) 



(18) 



with Ua and rif, the refractive indices of the signal and 
idler, and $(wfe) the phase matching function that takes 
into account the spatial matching of the modes in the 
crystal As usual, the frequency Vtk = [ujp — 

4^^]/2 has been introduced in the Hamiltonian 



to appropriately rescale the ground state energy. Notice 
that, in contrast to the case described in Sec. IHBI the 
Hamiltonian (|17l) couples nondegenerate modes whose 
frequencies sum up to the pump frequency ujp. Canonical 
transformations analogous to (|15|l allow us to rewrite the 
Hamiltonian in the free field form 



H ^ 



\h{ujk - ^k)AlAk + h{ujp -UJk- nk)BlBk (19) 



With this Hamiltonian, the partition function is given by 



InZ(A) = ^In 

fe 

+ El- 



1 



1 — g-A?i(wfc-nfc) 
1 



(20) 



where the two contributions are due to the signal and 
idler modes respectively and the sum over k is performed 
on all the frequencies up to ujp. For ease of calculation, 
we will assume that the coupling is acting only over 
a frequency band (^LUp (C < 1) centered around LOp/2, 
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where it assumes the constant value ^. This choice gives a 
rough approximation of the crystal phase matching func- 
tion $(a;fc) that prevents the coupling of signal and idle r 
photons when their frequencies are too mismatched 10]. 
In the high energy regime E ^ hSui, the sums in Eq. H20() 
can be replaced by frequency integrals, obtaining 



2 /•'-'p(i-C)/2 
InZ(A) = -- / duj \n 
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-Xhuj 



(21) 



process (amplification plus input energy), no gain is ob- 
tained. The correct attitude is, of course, to consider the 
signals alone. 
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where = [wp- y^cj^ - 4^2] /2. The integrals in Eq. if^ 

have no simple analytical solution, but we can give a 
perturbative expansion in the low-interaction regime i.e. 
e = 4^2 /cjp <C 1. In this limit, the result is derived in 
Add. 1X1 and is given by 



8lo 



Co 



f E5uj 



e ci 



c ) 



(22) 



where the functions cq and ci are plotted in Fig. [3 The 
zeroth order term cq in (|22() gives the capacity of two 
broad-band non-interacting modes with cutoff frequency 



in the limit of infinite bandwidth (w. 



this 



function reaches the asymptotic behavior C 

In 2 



00), 

Casym — 

of Eq. @ for 

To first 



y/E/ (3Mw) that corresponds to C^b 
two noninteracting wide-band bosonic systems 
order in e, the increase in capacity due to the interaction 
^ is given by the value of ci , which is a positive quantity 
(see App. 0. It can be shown that ci ^ in the limit 
Wp — ^ 00, so that in the infinite bandwidth regime no im- 
provement is obtained from the interaction. In fact, an 
infinite continuous spectrum is invariant under the trans- 



formation UJk 



As in the non-interacting 
broad-band case of Eqs. © and pOfl . the transmission 
time r of the signal can be estimated as 2tt/Slo. The 
validity of this assumption (which implies negligible dis- 
persion) rests on the fact that the Hamiltonian H17|l is 
valid to first order in the interaction term S,k and the 
small dispersion which derives from such term does not 
play any role in the capacity ^| . 

In this calculation E is the energy devoted to the sig- 
nal and idler modes: it does not include the energy spent 
to pump the crystal. This last quantity is proportional 
to l^pp (apart from corrections of order e which derive 
from the coupling). In the limit of undepleted pump- 
ing we have considered, this is a very large contribution 
which overshadows E and that is not directly used to 
encode information. In this respect, the system is very 
inefficient in using the total available energy. However, 
this is not the correct attitude to evaluate such a system: 
our process is analogous to the amplification of signals 
where, if one considers the total energy employed in the 



FIG. 2: a) Capacity function 00(7) (continuous line) and its 
correction £1(7, C^) (dashed line) of Eq. I22II with fractional 
coupling bandwidth C = -5 and 7 = ESui / {2hu)p). The capac- 
ity increase accomplished by the nonlinearity is evident from 
the positivity of the term ci . b) Comparison between the ca- 
pacity C of Eq. 1221 with e = =0.1 (continuous line) 
and the asymptotic two mode rate Casym = ^^\/27/3, 
obtained using an infinite frequency band (dotted line). 



D. Swapping Hamiltonian 

Consider N modes that are pairwise coupled through 
the Hamiltonian 

N 

H = nwa]aj + ^ Ajjva]ajv 

= fta^' • (wl-fA) - a, (23) 

where a is the column vector containing the annihilation 
operators aj of the N modes and A is an iV x iV sym- 
metric real matrix with null diagonal. This Hamiltonian 
describes TV modes aj whose photons have a probabil- 
ity amplitude Kjji to be swapped into the mode aji. By 
performing a canonical transformation on all the mode 
operators, it is possible to rewrite Eq. in the free 
field form 



N 



(24) 



where the A^s are the N eigenvalues of A. Two condi- 
tions must be satisfied: the positivity of H requires that 
Aj ^ Lo for all j, and, since the diagonalization of A must 
preserve its trace, Y^- Xj = 0. The capacity of this sys- 
tem can now be easily computed as P, |3, 



C = max 

ei,e2,---,eiv 



N 



h{uj + \j ) 



(25) 
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where the maximum must be evaluated under the re- 
quirement that Cj = E. In the simple case of = 2 
(where Ai = — A2 = S,), the maximization (|25|l can be eas- 
ily performed numerically: the increase in capacity over 
the two mode non interacting case is presented in Fig.|Sl 
In this case, in the strong coupling regime (|^| — > uj) the 
capacity diverges as log2[-E/ {h{uj— |C|))]: this corresponds 
to employing for the information storage only the lowest 
frequency mode among Ai and A2. 




FIG. 3: Capacity increase of the swapping Hamiltonian as a 
function of the energy parameter E/hw and of the coupling 
ratio given by plotting AC = C — 2g{E/ {2huj)) (the first 
term being the capacity C of Eq. I25II and the second the 
capacity of two non interacting modes that share an energy 
E). 



complicated, we show an increase in capacity by consider- 
ing the following choice of the coupling constants (which 
may well not be the optimal one) : 

Xjk ^ -ujkr for j = l,---,7V-l (28) 
XNk - LUk{N~l)r, (29) 

where ^ r 1. The choices (I2HI and (|2HJl auto- 
matically guarantee both the positivity condition on the 
Hamiltonian and the trace preservation condition on the 
A^'') matrices. With this choice, the Hamiltonian 1)2 7|l 
becomes 

Af-l 

H = (l-r)^^;i^feylt^yl,fe 

k 3=1 

+ [l + {N - I) r]Y,^kA\^ANk . (30) 

fe 

Essentially we have chosen a coupling that contracts by a 
factor 1— r the frequencies of the first A^— 1 normal- modes 
and stretches by a factor [1 -I- (A^ — 1) r] the frequency of 
the A^th one. Choosing r — > 1, we can increase the capac- 
ity ad libitum over the case of A^ non-interacting systems. 
In fact, in this limit a straightforward application of the 
wideband calculation of Eq. ©, gives C ~ ^^Cwb- 
This result must be compared with the case in which 
the A" wideband modes are independent where the ca- 
pacity scales as VlVCwb- Clearly an arbitrary increase in 
capacity is gained as r — > 1. 



E. Broadband swapping Hamiltonian 

A generalization of the Hamiltonian H23|) given in the 
previous subsection can be obtained by considering A^ 
parallel broadband modes in which the coupling joins all 
the modes with the same frequency, namely, 

N 

k ] = 1 

+ EEa5-?4«^'^' (26) 

k j^^j' 

where ajk is the annihilation operator of the jth system 
with frequency LUk- This Hamiltonian can be taken to a 
free field form using the procedure of the previous sub- 
section, so that 

N 

H^Y^T. ^(^fc + ^Jk)A],A,k , (27) 

k j=l 

where Xjk is the jth eigenvalue of the matrix A*^*^). Ex- 
pressed in terms of the normal modes Ajk, the Hamil- 
tonian (|27|l describes A^ independent wideband modes. 
Since in this case the maximization of the form H25|) is 



III. GENERAL CONSIDERATIONS 

In the previous sections we have derived the capacity of 
various types of nonlinear bosonic systems. The common 
feature of all these results is that the nonlinearities were 
used to reshape the spectrum by compressing it to lower 
frequencies, where it is energetically cheaper to encode 
information. A couple of remarks are in order. First of 
all, in our calculations the mean energy E represents the 
energy of the modes in the system, whereas customarily 
one considers the input energy. Our choice is motivated 
by the fact that nonlinear systems dispense energy to the 
information bearing modes, so that the input energy does 
not necessarily coincide with the amount of energy that 
the medium where information is encoded needs to sus- 
tain. This last quantity is a practically relevant one for 
the cases in which the degrees of freedom used to encode 
information cannot handle high energies (as for example 
optical fibers ^3)- Notice that, in some of the cases we 
have studied (the example of Sec. Ill Cl in particular), the 
nonlinearity is achieved by supplying the system with an 
external energy source in the form of an intense coherent 
beam. If one were to take into account also this contribu- 
tion in the energy balance, then no capacity enhancement 
would be evident, since the pumping energy is used only 
indirectly to store the information. However, it is not un- 
warranted to exclude the pump from the energy balance, 



6 



since it is only used to set up the required Hamiltonian 
and not directly employed in the information processing. 
Finally, our analysis is limited to the noiseless case. In 
the presence of noise in place of the von Neumann en- 
tropy in Eq. one would have to consider the Holevo 
information This is a highly demanding problem 

because of the yet unknown additivity properties of this 
quantity. At least in the case of linear bosonic systems, 
the capacity in the presence of noise was studied in |1.[T3|. 

This work was funded by the ARDA, NRO, NSF, and 
by ARO under a MURI program. 



band modes (signal and idler) with cut-off frequency Wp. 
Replacing Eq. l|Aip into the energy constraint we can 
find the value of the Lagrange multiplier A, contained in 
the parameter (3, by solving the equation 



dfoW , a/i(/3,c) 



8/3 



d(3 



(A4) 



where 7 = E5uj/{2fujj^) is a dimensionless quantity. By 
expanding the solution (3 for small e as /3 = /3o + f-Pi, it 
follows that 



APPENDIX A 

In this appendix we derive formula (|22|l for the broad- 
band parametric down-conversion Hamiltonian. In the 
low-interaction regime e <C 1, Eq. H21|l becomes 

lnZ(A)-^[/o(/3)+e/i(/3,0+0(e')] , (Al) 



where [3 = Xhujp and 
/o(/3) ^ 
/i(/3,C) ^ 
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1 dx In 


'0 


1 - e-^_ 



In 



1 _ g-/3(l+C)/2 



1 - 6-/3(1-0/2 



(A2) 
(A3) 



Notice that the zeroth order term /o in the expansion 
(|A1|) corresponds to the partition function of two broad- 



dfoWo) 



-7 



Pi 



dp 

dfi{Po,0 /d^MPo) 



dp 



dp^ 



(A5) 
(A6) 



These two equations can be numerically solved for any 
value of 7. Replacing the solution in ljAl|) and using 
the capacity formula © we can evaluate the parametric 
down-conversion capacity as reported in Eq. H22|l where 



co(7) = [/3o7 + /o(/3o)]/ln2 
ci(7,C) = /i(/3o,C)/ln2. 



(A7) 
(A8) 



Both these functions depend on the system energy only 
through the quantity 7. 
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